Heat Transport through Rough Channels 
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We investigate the two-dimensional transport of heat through viscous flow between two parallel 
rough interfaces with a given fractal geometry. The flow and heat transport equations are solved 
through direct numerical simulations, and for different conduct ion- convection conditions. Compared 
with the behavior of a channel with smooth interfaces, the results for the rough channel at low and 
moderate values of the Peclet number indicate that the effect of roughness is almost negligible on the 
efficiency of the heat transport system. This is explained here in terms of the Makarov's theorem, 
using the notion of active zone in Laplacian transport. At sufficiently high Peclet numbers, where 
convection becomes the dominant mechanism of heat transport, the role of the interface roughness is 
to generally increase both the heat flixx across the wall as well as the active length of heat exchange, 
when compared with the smooth channel. Finally, we show that this last behavior is closely related 
with the presence of recirculation zones in the reentrant regions of the fractal geometry. 

PACS numbers: PACS numbers: 47.27.Te, 44.15. +a, 47.53. -|-n 



I. INTRODUCTION 



The study of transport phenomena in irregular media 
is fully recognized nowadays as an important research 
subject with many technological and industrial applica- 
tions. In particular, the role of the morphology of an 
active interface on the efficiency of heat and mass trans- 
port devices has also been a theme of great interest for 
the scientific community. For example, the development 
of a modern approach for the project of heat exchang- 
ers must necessarily take the detailed description of the 
surface geometry as a possibility for the proper design 
of these devices. It is well-known in this field that the 
inclusion of extended surfaces or fins to the surface of 
the heat exchanger can substantially increase the rate of 
heat transfer in the system and, therefore, the efficiency 
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of the equipment 

For all practical purposes, the restrictions imposed to 
the geometrical features of an exchange interface should 
be the result of an intimate balance between its complex- 
ity and cost from one side, and the performance of the 
equipment from the other. In any case, an important 
issue in the design of heat exchangers is certainly the ex- 
tent to which the heat transfer capacity of the system is 
dictated by the available surface area. For instance, un- 
der diffusion-limited conditions, the access of heat to an 
irregular surface can be highly nonuniform due to screen- 
ing effects. In this situation, the regions corresponding 
to fins or extended protrusions will display a higher ac- 
tivity as opposed to the deep parts of fjords which are 
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more difficult to access. As a result, under strong dif- 
fusive constraints, the efficiency of the surface for heat 
exchange can be substantially different from the one ex- 
pected, considering its total area, the intrinsic conductive 
properties of the material, and the transport properties 
of the heating (or cooling) fluid. 

An analogous problem can be found in the context of 
electrochemistry. The nonuniform activity during the lin- 
ear transport through irregular interfaces of electrodes 
has been explained in terms of a screening effect which 
is entirely similar to the one observed in heat transfer. 
The extensive research developed on this subject in re- 
cent years H, ^ 1^ has been mainly devoted to the in- 
troduction, calculation and application of the concept of 
active zone in the Laplacian transport to and across irreg- 
ular interfaces. For example, through the coarse-graining 
method proposed in |3j| it is possible to compute the flux 
through electrode surfaces from their geometry alone, 
without solving the general Laplace problem. In a sub- 
sequent study 6], it has been shown that this technique 
provides consistent predictions for the activity of irregu- 
lar catalyst surfaces. More recently, the activity of irregu- 
lar absorbing interfaces operating under diffusion-limited 
conditions has been investigated through nonequilibrium 
molecular dynamics 0. The simulation results indicate 
that the extent of the interface that is significantly ac- 
tive is rather sensitive to the governing mechanism of 
transport. Precisely, the length of the active zone de- 
creases continuously with density from the Knudsen to 
the molecular diffusion regime. 

The study of laminar-forced-convective heat transfer 
through tubes or between parallel plates with smooth 
walls is well known in the literature as the Graetz prob- 
lem Q . Graetz solved the problem analytically with the 
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assumptions of steady, irrotational and incompressible 
flow, constant fluid properties, fully developed velocity 
profiles, and in the absence of energy dissipation effects. 
Much effort has been dedicated in the past years to im- 
prove the accuracy of the Graetz solution and generalize 
the Graetz problem to include other types of boundary 
conditions [1 113 IH 111 111 El El ■ The development 
and solution of more realistic models for heat transport 
in internal flow configurations and under forced convec- 
tion conditions has also been a research theme of great 
scientific interest Few studies, however, have 

been dedicated to the investigation of an important ex- 
tension of the Graetz problem, namely, the case in which 
the flow and heat transfer processes are confined to walls 
of irregular geometry. 

In the present work we apply the concept of active 
zone to investigate the steady-state transport of heat in 
a fluid flowing through a two-dimensional channel whose 
walls are irregular interfaces. We investigate through di- 
rect computational simulations the effect of the interface 
geometry on the heat exchange efRciency of the system 
for different diffusive (conductive)-convective conditions. 
Compared to the behavior of a smooth channel under 
conditions where the convective mechanism of transport 
is relevant, our results show that the activity of the in- 
terface can be significantly underestimated if geometrical 
details are not adequately considered. 

This paper is organized as follows. In Section 2, we 
introduce the general description of the physical system 
and study the flow through the irregular geometry. The 
transport of heat through the flowing fluid in the rough 
channel is then numerically treated in Section 3. Finally, 
general conclusions and perspectives are drawn in Sec- 
tion 4. 



II. VISCOUS FLOW BETWEEN ROUGH 
INTERFACES 



As shown in Fig. 1, our physical system is a two- 
dimensional channel of length L whose delimiting walls 
are identical interfaces with arbitrary geometry separated 
by a distance 2h. By definition, the two interfaces are 
made by the connection of several wall subsets. In the 
particular roughness model adopted here, we consider 
that these subsets are identical pre-fractal curves with 
the geometry of a square Koch tree [T^l . 

The mathematical description for the detailed fluid 
mechanics in the channel is based on the assumptions 
that we have a continuum, Newtonian and incompress- 
ible fluid flowing under steady state conditions. Thus, 
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Here the independent variables x and y denote the posi- 
tion in the channel, u and v are the components of the 
velocity vector in the x and y directions, respectively, and 
p is the local pressure in the system. The relevant phys- 
ical properties of the flowing fluid are the density p and 
the viscosity p. In our simulations, we consider non-slip 
boundary conditions at the entire solid-fluid interface. 
In addition, the changes in velocity rates are assumed to 
be zero at the exit x — L (gradientless boundary condi- 
tions), whereas a parabolic velocity proflle is imposed at 
the inlet of the channel. 



u{Q,y) = \v[l-{l) 



(4) 



The Reynolds number is defined here as Re = pVh/p, 
where V is the average velocity at the inlet. For simplic- 
ity, we restrict our study to the case where the Reynolds 
number is sufficiently low (Re w 1) to ensure a laminar 
viscous regime for fluid flow. Finally, we also assume that 
the fluid density and viscosity are constant properties of 
the flow throughout the channel. As a consequence, the 
flow field is decoupled from temperature and can be in- 
dependently calculated. 

The numerical solution of Eqs. Q-lEIl for the veloc- 
ity and pressure fields in the rough channel is obtained 
through discretization by means of the control volume 
finite-difference technique For the complex geom- 

etry involved, we solve this problem using a structured 
mesh based on quadrangular grid elements. For exam- 
ple, in the case of the channel with walls subsets that are 
3-generation Koch curves, a total of 79,200 cells gener- 
ates satisfactory results when compared with numerical 
meshes of small resolution. The integral form of the gov- 
erning equations is then considered at each cell of the nu- 
merical grid to produce a set of algebraic equations which 
are pseudo-linearized and solved using the SIMPLE al- 
gorithm I23. The criteria for convergence used in the 
simulations is defined in terms of residuals, i.e., a mea- 
sure of the degree to which the conservation equations 
are satisfied throughout the flow field. In all simulations 
we performed, convergence is considered to be achieved 
only when each of the residuals fall below 10~^. 

In Fig. 2 we show the velocity profiles corresponding to 
two different sections of the rough channel along its lon- 
gitudinal direction. Although very small when compared 
to the velocities at the center, the velocities of the fluid 
cells constituting the reentrant zones are still finite. In- 
deed, as depicted in Fig. 3, a close-up of these roughness 
details shows fluid layers in the form of consecutive eddies 
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whose intensities fall off in geometric progression |21ll22| . 
Although much less intense than the mainstream flow, 
these recirculating structures are located deeper in the 
system, and therefore experience closer the landscape of 
the solid- fluid interface. More precisely, viscous momen- 
tum is transmitted laterally from the mainstream flow 
and across successive laminae of fluid to induce vortices 
inside the fractal cavity. These vortices will then gen- 
erate other vortices of smaller sizes and intensities. In 
the next section, we show how this flow structure affects 
the overall transport of heat as well as the distribution 
of fluxes among the active wall elements in the rough 
channel. 



III. HEAT TRANSPORT 

Once the velocity and pressure fields are obtained for 
the flow in the rough channel, we proceed with the mathe- 
matical modeling and numerical computation of the tem- 
perature scalar field. As already mentioned, we assume 
that the flow is independent of the temperature. In prac- 
tice, we simply "freeze" the velocity field solution previ- 
ously calculated and use it to compute the steady-state 
temperature field T(x,y) solving the two-dimensional 
diffusion-convection equation. 



dT dT 

+ — + — h (5) 



where a = k/pCp is the thermal diffusivity, k is the ther- 
mal conductivity, and Cp is the specific heat of the fluid. 
For boundary conditions, we consider that the inlet is 
subjected to a constant temperature Tq and that the tem- 
perature Tuj is imposed on both wall interfaces. In partic- 
ular, we study the case in which Tq > T^,, i.e., the heat is 
transported from the bulk of the fluid to the walls of the 
rough channel. The relative influence of the conductive 
(diffusive) and convective mechanisms of heat transport 
on the behavior of the system is quantified here in terms 
of the Peclet number, defined as 



Pe = Vh^/aL 



(6) 



Similarly to the procedure already utilized for the calcu- 
lation of the flow, the temperature T{x, y) is also ob- 
tained through discretization of Eq. |(SJl by means of 
finite-differences, with an upwind scheme to avoid nu- 
merical instabilities due to the presence of convection 
|2Clj |. Using this numerical technique, simulations have 
been performed for three types of channels correspond- 
ing to the first three pre-fractal generations of the square 
Koch tree. 

In Fig. 4 we show contour plots of the temperature 
field in a 3-generation channel for four different diffusion- 
convection conditions. At very low Pe (Fig. 4a), conduc- 
tion is the dominant mechanism with the transport of 
heat being practically limited to the entrance region of 
the channel. By increasing the value of Pe (Figs. 4b-d), 



the contribution of convection to heat transfer through 
the center of the channel becomes relevant. As shown 
in Fig. 4c, although the inlet temperature can reach the 
exit at Pe « 200, the temperature at the reentrant zones 
of the wall remains virtually unaffected by convection. 
Only at very high values of Pe, when the vortices in- 
side these irregularities can effectively transport heat, the 
high temperature front can penetrate deeper in the sys- 
tem (Fig. 4e). 

The role of the wall roughness on the transport of heat 
in the channel can be quantified in terms of the normal- 
ized heat flux 
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where q^ is the total heat flux crossing the fluid-solid 
interface from the bulk. 
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The reference value go corresponds to the heat flux cross- 
ing a smooth interface of size L at temperature Tu, that is 
part of a system where heat is transported by conduction 
through a static fluid from a line at a fixed temperature 
To, 
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(9) 



The dependence of the flux on the Pe number is shown 
in Fig. 5 for the flrst three generations of the selected pre- 
fractal geometry. The behavior of a smooth channel of 
size L and subjected to the same boundary conditions 
is also shown for comparison. In this simpler case, the 
system displays a sharp increase of the flux with Pe, up to 
a point (Pe « 10) where the following power-law regime 
persists for more than 5 orders of magnitude: 



Pel/3 



(10) 



This behavior has been observed in previous theoretical 
studies In the case of the power-law, it has been 

demonstrated that the analytical solution of the classical 
Graetz problem, where the mechanism of axial conduc- 
tion is disregarded, converges asymptotically to this scal- 
ing behavior |l3j |. Alternatively, as we show in the Ap- 
pendix, this typical power-law regime of heat transport, 
known to be valid for moderate and high Pe values, can 
be readily recovered by means of a global heat balance 
on the system combined with a simple scaling argument. 

Remarkably, the results shown in Fig. 5 indicate that 
different rough channels display approximately the same 
behavior for low and moderate Pe numbers (Pe < 200). 
Moreover, their performances are also similar to the be- 
havior of the smooth channel for the same range of Pe 
values. Under these conditions, the conduction of heat 
either predominates over convection or still represents 
a significant contribution to the overall fiux across the 
walls. 
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For Peclet values in the range 10 < Pe < 200, each 
pre-fractal subset of the interface with the surrounding 
fluid can be approximately taken as a purely diffusive 
cell subjected to a constant temperature at the wall, 
and a temperature Tq at a line in the bulk of the fluid 
separated by a given distance from the center of the chan- 
nel. Considering steady-state conditions, the system can 
then be treated as a two-dimensional Laplacian transport 
unit subjected to Dirichlet's boundary condition. This is 
a typical situation where the the theorem of Makarov 
holds p^ . This theorem states that the information di- 
mension of the harmonic measure on a singly connected 
interface in d = 2 is exactly equal to 1, where d represents 
dimension. Makarov's theorem has a simple but impor- 
tant consequence: the screening effect due to geometrical 
irregularities of the interface can be accessed in terms of 
the ratio S = Lp/ L 3], where Lp is the perimeter of the 
interface. Because the active length of the interface La is 
of the order of the system size, La ~ then La ~ Lp/ S 
and the factor S can be considered as the "screening fac- 
tor" of the Dirichlet-Laplacian field. When applied to 
our physical system, the practical meaning of this exact 
result is that, irrespective of the surface morphology, the 
length of the zone within the heat exchange wall which 
receives most of the heat flux should be of the order of the 
system size L, instead of the perimeter Lp. This explains 
why in Fig. 5 the heat flux at the interface is nearly in- 
sensitive to differences in the roughness geometry for low 
and moderate values of Pe. 

For Pe > 200, the predominance of convection over 
conduction is extended to the reentrant zones of the pre- 
fractal units. This crossover value results from the fact 
that the magnitudes of the velocity vectors constituting 
the largest vortices in the rough region are approximately 
two orders of magnitude smaller than the main flow ve- 
locity V. At high Pe values, heat can then be convected 
through these recirculating structures of the flow towards 
the less accessible wall elements of the interface. Under 
such circumstances, the theorem of Makarov is not valid 
anymore, and the flux (p across the wall becomes depen- 
dent on the geometrical details of the interface. As shown 
in Fig. 5, the higher the generation of the pre-fractal in- 
terface, larger is the heat flux at high Peclet conditions. 

We now show that the previous analysis based on the 
overall heat flux </) is consistent with the behavior of a 
more detailed measure of transport. In particular, the 
flux heterogeneity at the interface can be quantified in 
terms of an active length defined as 

Lp 

La = l/Y.(fi {l<La<Lp), (11) 
i=l 

where the sum is over the total number of interface el- 
ements Lp, (pi = Qi/'^Qj, and qi is the thermal flux at 
the wall element i. From the definition (|ll|l . La — Lp 
indicates a limiting state of equal partition of normalized 
fluxes ~ 1/Lp, Vi) whereas Lq = 1 should correspond 
to the case of a fully "localized" flux distribution. The 



active length therefore provides an useful index to quan- 
tify the interplay between the Laplacian phenomenology 
and the complex geometry of the absorbing interface at 
the local scale. As expected, the results shown in Fig. 6 
indicate that La generally increases with Pe. It is im- 
portant to note, however, that the activities of smooth 
and rough interfaces behave approximately in the same 
way for low values of Pe. This is also consistent with 
Makarov's theorem under conditions where the conduc- 
tion mechanism controls the transport of heat inside the 
pre-fractal roughness. 

For high Peclet values. La becomes strongly sensitive 
to the geometrical details of the interface. Again, due to 
the relevant contribution of convection to the transport 
at the smaller scales, the flux distribution among the in- 
terface units becomes more uniform. The higher the pre- 
fractal generation used to create the roughness geometry, 
higher is the available number of wall elements for heat 
exchange. As a consequence, a substantial increase in the 
active length can be observed for a fixed value of Pe. 



IV. CONCLUSIONS 

The development of a coherent scientific approach to 
the project of heat and mass transport devices involving 
new concepts from mathematics and physics still repre- 
sents an important research and technological challenge 
to chemical and mechanical engineers. In the present 
work, we have studied the conditions under which the 
roughness geometry of a channel with a fiowing fiuid 
should play a relevant role on the transport of heat. For 
low Peclet numbers, when the mechanism of heat con- 
duction dominates, it is shown that the overall flux of 
heat across a pre-fractal interface subjected to Dirichlct 
boundary conditions has little sensitivity to details and 
irregularities of the heat absorption zone. Under similar 
conditions, the length of the active zone computed for 
a smooth channel is very close to the measure obtained 
for a rough system. This behavior is then interpreted in 
terms of a screening phenomenon that privileges the ab- 
sorbing activity of the geometrically exposed sites of the 
interface. Moreover, such a "localization" of the heat flux 
can be formally explained in terms of Makarov's theorem 
for Laplacian transport in two-dimensional irregular ge- 
ometries. 

The situation is entirely different at high Peclet num- 
bers. In this case, our simulations reveal that both the 
heat flux and the active length start to depend noticeably 
on the morphological details of the interface. In this sit- 
uation, not only the total heat flux leaving the system 
increases faster for higher generations of the pre-fractal 
interface, but the distribution of fluxes at the wall units 
also becomes gradually less "localized" {La increases) as 
Pe increases. This transition refiects the onset of con- 
vective effects on the heat transport in the proximity of 
the rough interface. A closer look at the system indi- 
cates that the significant changes observed in La above 
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the transition point are induced by the presence of recir- 
culation zones of flow (vortices) in these regions. 

As a future work, we will extend the approach intro- 
duced here to simulate the effect of different geometries 
on the transport efficiency of the system as well as other 
types of "absorption" mechanisms (e.g., finite-rate heat 
transfer at the walls) limited by diffusion transport. Fi- 
nally, we also believe that the results presented in this 
study provide useful information to the understanding of 
other physico-chemical systems with relevance to science 
and technology that include heat, mass or momentum 
transport. The modeling strategy devised here should 
also be helpful as a design tool to choose a suitable in- 
terface geometry for a given heat exchange system. 



V. ACKNOWLEDGMENTS 

We thank Bernard Sapoval and Marcel Filoche for 
many helpful discussions. We acknowledge the Brazil- 
ian agencies CNPq, CAPES and FUNCAP for financial 
support. 



VI. APPENDIX: THE SMOOTH CHANNEL 

Here we revisit the classical problem of heat trans- 
port through a fluid flowing in a channel whose walls are 
smooth interfaces. The system setup is shown in Fig. 7. 
As in the case of the rough channel, we consider a New- 
tonian and incompressible fluid with constant properties 
that is flowing under viscous and steady-state conditions. 
Assuming non-slip boundary conditions at the walls, the 
velocity along the channel follows the classical parabolic 
profile 



(12) 



where i] = y/h, V = Ap h/{3^L), and Ap is the pres- 
sure drop between the entrance and the exit of the chan- 
nel. Under these simple flow conditions, the fluid is sud- 
denly subjected to the temperature Tq at the channel in- 
let, while a constant temperature T^, is imposed at both 
walls. At steady-state, and assuming that the axial heat 
conduction mechanism is negligible, this corresponds to 
the so-called Graetz problem, where the temperature field 
in the channel obeys the following differential equation: 



with boundary conditions 

T{0,y)^To and T{x,h)^T^ 



(13) 



(14) 



An analytical solution for the linear partial differential 
equation H13|l subjected to H14|l can be obtained through 
separation of variables and a proper expansion in terms of 



the so-called Graetz functions From this solution, 

it is possible to show that the overall heat fiux through 
the absorbing walls, (f), follows the power law behavior, 
(j) ^ Pe^/"^, that is valid for moderate and high Pe values. 

The objective here is to show that this power-law 
regime of heat transport can be readily recovered by 
means of a global heat balance combined with a sim- 
ple scaling argument. First, in the absence of sources or 
sinks of heat, the total heat crossing the walls should be 
equal to the difference 



lout 



(15) 



where qin and Qout are the fiuxes at the entrance and at 
the exit sections of the channel, respectively. At suffi- 
ciently high values of the Peclet number, the convective 
part of the flux dominates the heat transfer process. As 
shown in Fig. 8, the temperature profile at the exit dis- 
plays a central plateau in the range —y* < y < y* , where 
the temperature is equal to the entrance temperature Tq, 
and then decays monotonously towards the wall temper- 
ature Tuj, in the region y* < \y\ < h. The value of 
?7*(= y* /h) can be found from the definition of the local 
Peclet number 



Per, = 



TC 



L/u{ri) 



(16) 



where td = 5^ /a is the characteristic conduction time, 
Tc = L/u is the local convective time, and 5 = h{l — rf). 
If we now argue that 77* should correspond to the position 
where Pe,, « 1, we can write from Eqs. H12|) and (|16|) that 



(l-r;*)^[l-(r;*)1 = l. 



2 aL 

which, from the definition leads to 

2 



(l-77*)^(l+77*) 



3Pe 



For high values of Pe, rj* is very close to 1, such that 
(1 — ri*Y{l + -q*) K, 2(1 - 77*)^, resulting in 



ri w 1 



1 



(3Pe)i/3 



(17) 



Going back to Eq. H15|) . the heat flux q^n entering the 
system is given by 

qin = 2/9CpTo / u{y) dy . 
Jo 

The flux gout leaving the channel can be expressed as the 
sum 



qout = qi+q2 



(18) 



where the first term corresponds to the heat escaping 
through the plateau region of temperature Tq and size y* 



(19) 



qi = 2pCpTo / u{y) dy , 
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while the second term represents the flux through the 
remaining part of the channel cross-section with variable 
temperature 



(?2 = 2pcp / u{y)T{y) dy 

Jy- 



(20) 



As a first approximation, if we assume that the temper- 
ature in this region decreases linearly with y, 



T{y) 



y* ~h 

Eq. H20() can then be written as 

92 = 2/9CpTo / u{y)dy 
J y 



{y-y*) + To 



+ 2pcp / u{y) 



ir y 

and the total flux becomes 



{y-y*)dy , 









['-©1 







{y - y*)dy . 



Rewriting this expression in terms of the dimensionless 
variable ry, we obtain 



^ V J rj* 



From this equation and the deflnitions © and Q, the 
normalized flux can be expressed as 



and using the approximation H17|l . we finally get the scal- 
ing relation 



(3Pe) 



1/3 
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(21) 
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FIG. 1: Schematic diagram of the rough channel under study. 
The wahs are a sequence of irregular interfaces with the geom- 
etry of squaxe Koch trees. The fluid flows steadily from left 
to right at Re m 1. A constant temperature To is maintained 
at the inlet, while the temperature on both walls is T^. We 
study the steady-state transport of heat for the case in which 
To > T^, i.e., the heat goes from the bulk of the fluid to the 
walls of the channel. 



FIG. 2: Longitudinal velocity profiles at two distinct cross- 
sections of the rough channel. Although small when compared 
to the mainstream flow, the velocities inside the reentrant 
zones of the pre-fractal geometry are not equal to zero, as can 
be observed in the profile on the left. At the cross-section on 
the right, the velocity profile is typically parabolic. 
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FIG. 3: Vortices in the reentrant zones of the pre-fractal 

roughness. A closc-up view of a smaller roughness detail re- 
veals the sequence of eddies with intensities that fall off in 
geometric progression. 



FIG. 4: Contour plot of the temperature field for different 
values of the Peclet number: (a) Pe = 0.25, (b) Pe = 10, (c) 
Pe = 200, (d) Pe = IQr' , and (e) Pe = 10^. The temperature 
decreases from red (gray dark) to blue (gray light). 
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FIG. 5: Log-log plot showing the dependence of the normal- 
ized flux across the walls, 0, on the Pcclot number, Pe. The 
solid line represents the behavior of a smooth channel, while 
the dashed, dot-dashed and dotted lines correspond to the 
first, second and third pre-fractal generations, respectively. 
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FIG. 6: Logarithmic plot showing the variation of the length 
of the active zone, La, with the Peclet number, Pe. The solid 
line corresponds to the case of a smooth channel, while the 
dashed, dot-dashed and dotted lines represent the behavior 
of the first, second and third pre-fractal generations, respec- 
tively. 
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FIG. 7: Schematic diagram of the smooth channel showing 
the fully developed parabolic profile of the flow. As in the 
case of the rough channel, the inlet is at temperature To while 
the temperature on both walls is maintained at Tw We also 
consider the case in which To > T^. 




FIG. 8: Steady-state profiles of the normalized temperature 
at the exit of the smooth channel for three different values of 
the Peclet number. 



